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Abstract 

The explanation of existence and ubiquity of 1/f- noise is under considera- 
tion based on the idea that 1 /f-noise has no relation to long lasting processes 
but originates from the same mechanisms what are responsible for the loss of 
correlations with the past, shot noise and fast relaxation. According to this 
idea, in a system which produces incident events being indifferent to their 
totally happened amount the fluctuations in the production of events grow 
proportionally to its average value. Such a free memoryless behaviour re- 
sults in the 1/f long-living statistical correlations without actual long-living 
causality. 

The phenomenological theory of memoryless random flows of events 
and of related Brownian motion is presented which closely connects 1/f- 
noise and non-Gaussian long-range statistics, both expressed in terms of 
only short-range characteristic scales. 

The exact relations between 1/f- noise and equilibrium four-point cumu- 
lants in thermodynamical systems are analysed. The presence of long-living 
four-point correlations and flicker noise in the Kac's ring model is demon- 
strated. 

The problems of kinetical approaches to noise and irreversibility in 
Hamiltonian systems are touched, and the general idea is confirmed in the 
case of gas. 

It is argued that correct construction of gas kinetics in terms of Boltz- 
mannian collision operators needs in the ansatz whose meaning is conser- 
vation of particles and probabilities at the path from in-state to out-state 
inside the collision region. This reason results in the infinite set of kinetical 
equations which describe the evolution of many-particle probability distri- 
butions on the hypersurfaces corresponding to encounters and collisions of 
particles. These equations reduce to usual Boltzmann equation only in the 
spatially uniform case, but in general forbid molecular chaos because of 
spatial correllations of colliding particles. 

The formulated kinetics are applied to investigate the statistics of self- 
diffusion in equilibrium gas. The peculiar behaviour of the four-order cu- 
mulant of Brownian displacement of a gas particle is found being identical 
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to 1/f- fluctuations of diffusivity and mobility. This is the example of dy- 
namical system which produces no slow processes but produces 1/f-noise. 
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I. INTRODUCTION 

The low-frequency 1/f-noise (flicker noise) was discovered by Johnson in 1925 in elec- 
tronic lamps. At present 1/f-noise is known in giant variety of systems and processes 
in physics, chemistry, biology, medicine, cosmic and Earth sciences and in human activ- 
ity. Perhaply, thousands of experiments were performed during many years in order to 
comprehend this apparently mysterious phenomenon. Howevere, a conventional theory of 
1/f-noise is still lacked. 

It seems that the situation testifies the absence of true ideas. May be, that is why 
there are no new substancial reviews on 1/f-noise like ones [1,2] published in 1981. As 
before, usually 1/f-noise is believed be produced by some slow processes with a broad 
distribution of long relaxation times (life-times), and 1/f spectrum is thought as super- 
position of Lorentzians. But in most particular cases the hypothetical slow processes are 
not identified. 

The popular idea about long life-times is evoked by common opinion that any long- 
range statistical correlations associated with 1/f spectrum reflect some real long-living 
causal connections between the events which took place in the past and which could 
occur in the future. But as long ago as in 1950 N.S.Krylov [3] argued the necessity to 
make difference between statistical correlations and factual causality. 

In view of this distinction, the non-trivial explanation of 1/f-noise is possible [4-11] 
attributing its origin just to absence of long-living causality. More concretely, to processes 
which currently forget their own history and thus could not keep an eye on its character- 
istics. The flow of incident events produced by such a process have no certain "number of 
events per unit time", and just this circumstance implies 1/f-noise. Indeed, if a process 
remembers nothing about old events and thus about previous number of events per unit 
time, how could it know about the future one? If the past is of no importance, all the 
more the future is too. 

From this viewpoint, the long-range 1/f correlations are only the manifestation of a 
freedom of the random process. Correspondingly, the shorter is life-time of real causal 
memory the larger is magnitude of 1/f fluctuations [4-8], i.e. just the inverse rule takes 
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place, and only if the number of events per unit time is somehow anticipated by underlying 
dynamics one gets a process without 1/f- noise. 

For a first look, this picture is in contradiction with traditional hopes that absence 
of long-living memory should result in white noise, Marcovian behaviour and Poissonian 
flows of events. But in fact the latter is only very special variant of pure randomness. 
One has not to break foundations, merely the central limit theorem is not valid, and so 
another chapters of the probability theory [12] and theory of randomness in dynamical 
systems [13] should be involved. 

The main purpose of this maniscript is to review our own works on this subject because 
their English translations are hardly available. We do not try to review either the infinity 
of experimental data or recently suggested models of 1/f-noise (although some models may 
be in fact more close to our ideology than to conventional tendency to long life-times). 
Instead, we would like to show that if one avoids the limitations and approximations stan- 
dardly attracted when constructing models of noise then both the phenomenology [4-8,10] 
and approaches based on statistical mechanics [9,11] naturally predict 1/f contribution to 
whole noise without long life-times and slow processes. 

In particular, the correct procedure of derivation of gas kinetics from exact BBGKY 
equations demonstrates [9] that 1/f-noise exists even in gas in the Boltzmann-Grad limit, 
that is in the system where nothing like slow processes can exist. 

Unfortunately, up to now those who deal with rigorous statistical mechanics payed al- 
most no attention to 1/f-noise and likely had no suspicions that its nature could be as fun- 
damental as that of relaxation and irreversibility in Hamiltonian systems. In our opinion, 
1/f-noise is typical companion of relaxation and dissipation phenomena in time- reversible 
dynamics, and its sense is the absence of absolutely certain (well time-averagable) quan- 
titative rates of relaxation and irreversibility. 



Let X(t), (X) and SX(t) = (X(t) — (X))/ (X) be a random process, its mean value 
and its relative variation, respectively. The process will be called 1/f-noise if the power 
spectrum Ssx(f) of SX(t) looks approximately as 1/f at low frequencies and has no 
tendency to saturation at / — > . At high frequences, in practice Sgx(f) ~ * 5oo 7^ (the 
white noise level). The brackets () denote the ensemble averaging. Sometimes spectra like 
Ssx(f) f 1 are observed with 7 noticably different from unit, but we shall not discuss 
this unprincipal difference, except the meaning of formal non-stationarity at 7 > 1. 

The distinctive temporal property of 1/f-noise is the impossibility of precise measuring 
its mean value (X) with the help of time smoothing procedure. Indeed, the relative 
deviation 



II. THE PROBLEM OF 1/F-NOISE 



A. FORMAL DEFINITION OF 1/F-NOISE 




of the time- averaged value from (X) has the variance 
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(5x 2 (T)) » J [sin(7r/T)/7r/T] 2 ^(/)rf/ 



with /o being the inverse duration of observations. When T increases, first this integral 
decreases as {5x 2 (T)) pa Soo/T , in accordance with the large number law, but later turns 
into nearly constant. 

Formally, Sx(T) could either decrease in an extremely slow logarithmic way or even 
slowly grow if 7 > 1 . But practically any case results in existence of nonzero limit of 



y (5x 2 (T)) at large intervals T . Therefore, some characteristical uncertainty of the mean 
value of X(t) does exist, the so-called 1/f-floor or flicker floor. 



a) Let X(t) be the number of alpha-particle decays per unit time in a piece of radioac- 
tive matter (for instance, in Am with life-time of order of 1000 years). The experimentally 
discovered 1/f fluctuations of the rate of decays give the evidence that the decay proba- 
bility possesses relative uncertainty ~ 10~ 4 . 

b) Let X(t) be the electron drift velocity in a current-carrying slightly doped semicon- 
ductor. In accordance with the empirical Hooge's formula, S$x(f) ~ j , where a ~ 2-10~ 3 
, at frequences lower than tens kHz [1,2]. The corresponding relative uncertainty of drift 
velocity and mobility is of order of 0.1 , i.e. surprisingly large. 

c) The giant relative uncertainty, ~ 1 or even larger, characterizes 1/f-fluctuations of 
the inner friction and dissipation rate in excited quartz crystals [14,15]. In this case X(t) 
can denote the energy dissipated per unit time. 



i) Different solid, liquid and gaseous electric conductors [1,2]. X(t) is the charge 
transported per unit time, or resistance. Many of experiments are well described by 
Ssx(f) ~ ct/nf , where a = 10~ 7 — 10 1 for various media and n is the number of charge 
carriers in a sample under observation. 

ii) Highway traffic. X(t) is the number of cars passing per unit time. 

iii) Laser light Raman scattering [16]. X(t) is the number of photons scattered per 
unit time. 

iv) Chemical reactions in small volumes. X(t) is the number of reaction events per 
unit time. 

v) Definite type of diabetics. X(t) is the consumption of insuline per unit time. 

vi) Cosmic rays. X(t) is the current intensity of cosmic radiation. 

vii) Music. X(t) is the volume or dominating tone or other property of a sound track. 

viii) Earth's rotation. X(t) is the angular velocity of rotation. 

ix) DNA. X(t) is the number of definite complementary pairs per unit length of DNA 
sequence. In this case the distance along the sequence serves instead of time (for review 
and some original ideas see [17,18]). 

x) Thermal equilibrium voltage noise of a currentless disconnected conductor. X(t) is 
the instantaneous power spectrum S of the noise (that is the spectrum currently obtained 




B. EXAMPLES OF THE FLICKER FLOOR 



C. IT IS EVERYWHERE! MORE EXAMPLES 
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by filtering, squaring and smoothing over a finite time interval). Voss and Clarke [19] 
supposed that the resistance 1/f- noise in nonequlibrium current- carrying conductors (the 
example i) ) means the existence of similar 1/f- fluctuations of the power spectrum S of 
equilibrium white noise, in accordance with Nyquist formula if extended to slowly varying 
resistivity, and experimentally confirmed their hypothesis. 

D. 1/F-NOISE OF EQUILIBRIUM WHITE NOISE 
AND HIGHER-ORDER CUMULANTS 

The example x) is of special importance for us, because it demonstrates that 1/f-noise 
in disturbed thermodynamical systems can be theoretically investigated with the help of 
four-point equilibrium correlators. Indeed, the spectrum S of equilibrium noise u(t) is 
measured as a quadratic functional of u(t) , so the correlation function and spectrum of 
fluctuations of S can be expressed with four-point correlators of u(t) [20]. 

To be more precise, the information about 1/f-noise is hidden in those specific part 
of four-point correlators which can not be reduced to lower order correlations, i.e. in the 
four-point cumulants (semi-invariants) (^1,^2,^3,1x4) defined by 

(M1M2M3M4) = (U1U2) (U 3 U A ) + (uiU 3 ) (U2U4) + (U1U4) (u 2 u 3 ) + (u 1 ,u 2 ,u 3 ,u 4 ) 

where Ui = u(ti) and the equaliy (u(t)) = is supposed. In most practically interesting 
cases the equilibrium noise u(t) looks as white noise with a single characteristic relaxation 
time r c which determines its average spectrum, (S) = J (u(r)u(O)) dr m 2 (u 2 ) r c . The 
first three terms on the right-hand side reflect only this short-range relaxation and say 
nothing about low-frequency fluctuations of S described by the last cumulant term. 

If the equilibrium noise u(t) was purely Gaussian, this term would turn to zero. So the 
existence of 1/f-noise of white noise inevitably means that u(t) possesses non-Gaussian 
statistics. 

E. WHAT IS THE QUESTION? 

Though being everywhere the 1/f-noise phenomenon still stays without a commonly 
accepted explanation. Therefore the natural question arises: could it be explained in 
some general manner or not? 

At present, the most popular approaches to this issue are based on the idea that 1/f- 
noise is caused by some "slow processes" having large life-times [1,2]. For instance, in 
case of electric 1/f-noise in solids the slow processes could be associated with scattering of 
electrons by metastable movable structural defects or with their trapping by long-living 
structural inhomogenities. With no doubts, some part of low-frequency noise in solids is 
produced by such mechanisms and reflects real slow processes. It is interesting that formal 
summation of Lorentzians corresponding to a distribution of activation energies can re- 
produce factually observed connections between frequency and temperature dependencies 
of 1/f-noise in metals [1]. 

But, at the same time, there are no doubts that this aproach is unsufncient for complete 
explanation of 1/f-noise even in solids, not to mention the electric 1/f-noise in liquid 
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metals and other liquids, all the more, various 1/f-noises as in the above examples. May 
me, Tandon and Bilger [21] were right that the nice idea about summation of Lorentzians 
perhaply is a disaster for building a true 1/f-noise theory. 

III. PHENOMENOLOGY AND PHILOSOPHY 
OF 1/F-NOISE 

A. 1/F WITHOUT SLOWNESS 

The non-standard idea we are carrying from 1982 [4-11] is that 1/f-noise in general 
has no relation to actual (physical, chemical, etc.) slow processes, long-living memory 
and causal correlations, but, in opposite, manifests the absence of thats. 

To feel the idea let us consider the example v) . Assume that after taking a doze of 
insuline a diabetics completely restores the basic normal state of his organism and takes 
next doze only under the necessity, not under a time-table. By definition, the basic state 
does not depend on how frequently he taked the drag during last day, last week, etc., all 
the more during last year. The past is constantly forgotten. If it is so, then the organism 
has neither reasons nor ways to establish some definite " doze per unit time" . And one has 
no grounds to hope that the time avering will bring him some predictably certain result, 
just because of absence of long-living memory! The events are accumulated following the 
principle "let it be what occurs to be". The future is as unknown as the past. 

But such a free behaviour, without certain limit of the time smoothed doze, means 
the existence of 1 /f-fluctuations of "doze per unit time". Why 1/f , not other? If the 
smoothed doze had some definite limit but dependent on occasion, the low-frequency 
spectrum of the doze would acquire the delta-function 5(f) . In absence of the limit the 
spectrum should have the same dimensionality on the frequency scale, that is be nothing 
but 1/f [4]. 

By its nature, the 1/f-noise produced in such a way has no saturation at zero frequency, 
and one can not saturate it by a noisy perturbation. 

B. STATISTICAL LONG CORRELATIONS 
WITHOUT LONG-LIVING CAUSALITY 

The existence of the 1/f floor means that in a random flow of events the probable 
random deviation of the number of events from its ensemble-averaged value grows nearly as 
the average, the deviation is approximately proportional to time. The loss of information 
about old events gives the natural ground for this property : merely the system can not 
distinguish what part of happened events could be qualified as a proper part and what as a 
deviation from it. Any present result equally may serve as the starting point for further 
creation of events. 

The remarkable paradox is that this memoryless behavior can not be statistically de- 
scribed somehow except the language of long-living correlations! 
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All the events which belong to one and the same time piece seem mutually correlated, 
because all thats equally participate in the result. But there are no actual long-living 
causality beyond the corresponding statistical correlations [4-10]. 

Thus slow processes and slow statistical correlations are not one and the same thing. 
This important circumsance was underlined by N.Krylov [3] in 1950, in the frame of his 
critical revision of statistical thermodynamics. But he pointed out no examples. We 
supposed in 1982 that 1/f-noise is just the case [4]. 

C. SELF-FORGETTING FLOWS OF EVENTS 
AND CAUCHY STATISTICS 

Let us consider a stationary random process of accumulation of some quantity Q(t) 
, with instantaneous speed u(t) and AQ(t) = Q(t) — Q(0) = $Qu(t')dt' being the in- 
crement during time interval (0, t) , t > . For instance, AQ may be an amount of 
drag consumpted by diabetics. In other example AQ is spatial displacement of Brown- 
ian particle subjected to a constant force x [6]. In this case we take in mind so large t 
that the drift component of displacement much exceeds its diffusional component, that 
is Uot » \J1D§t (here Uq , -Do and \x = D$/T are the average drift velocity, diffusivity 
and mobility, defined by the relation Uq = (u(t)) = fix, with u(t) and T being velocity 
and temperature, respectively). One more example is the charge transport in electric 
junction, if applied voltage much exceeds T/e and so almost all electrons jump in one and 
the same direction. Then the number of transported electrons AQ(t) also can be treated 
as a random flow of events (random walk) with definitely (let positively) directed steps, 
AQ(t) > . 

Let To be the maximum of time scales characterizing causal correlations of u(t) , i.e. the 
life-time of dependence of the process on its history. We want to analyse the asymptotical 
statistics of AQ(t) and of the time-averaged speed (rate) of the walk U = AQ/t (doze 
per unit time, drift velocity, current, etc.) at t » r , when all the real memory effects 
should disappear [6,8]. 

Consider the amounts accumulated during two successive time intervals each longer 
than t , AQ(t) = Q(t) - Q(0) and AQ'(t) = Q(2t) - Q(t) , the summary amount 
AQ(2t) = AQ(t) + AQ'(t) and the corresponding time-smoothed rates U(t) = AQ(t)/t , 
U'(t) = AQ'(t)/t and U(2t) = AQ(2t)/2t = [U(t) + U'{t)}/2 . 

Because of t » r we believe that U and U' are independent random values. From 
the other hand, at t » r the result should be insensitive to t , because the system 
forgets how long ago the measurement had started and how large is t as compared with To 
. Consequently , U(2t) should coinside with U(t) in the statistical sense [6], i.e. 2U(2t) = 
U(i) + U'(t) is equivalent to 2U(t) —U(t) + U(t) . We come to conclusion that the sum 
U + U' of two independent random values behaves as the sum U + U of two completely 
dependent values] In other words, AQ(2t) does not statistically differ from 2AQ(t) . 

That is nothing but the distinctive feature of Cauchy statistics. It is easy to verify 
this in case of a random value z subjected to symmetrical Cauchy distribution whose 
probability density is W(z,w) = ^(w 2 + z 2 )^ 1 , with w being the width of distribution. 
Indeed, the convolution of W(z,w) with itself equals to W(z,2w) which is nothing but 
the probability density of the doubled variable 2z . 
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The same is true also with respect to maximally asymmetrical Cauchy distribution [12] 
whose characteristic function (CF), / exp(ikz)W(z, w)dz , looks as exp[— ikw ln(— ikcw)} 
, where the width parameter w determines the magnitude of random value, z oc w , and 
c ~ 1 . 

Obviously, if z is identified with AQ then the width w must be proportional to time, 
w oc t . Hence, the statistical equivalence of U (2t) and U (t) determines both the proba- 
bilistic and temporal properties of AQ(t) . 

Because Cauchy statistics makes no difference between dependent and independent 
contributions to a random walk, it is the naturally suitable statistics to describe correla- 
tions what have no actual causality beyond them. Its characteristical temporal property 
is that the magnitude of accumulated fluctuations is proportional to time, AQ(t) oc t 
(that is the central limit theorem fails). 

D. CAUCHY STATISTICS AND LONG CORRELATIONS 

Now we come to important point. The formal defect of the ideal Cauchy statistics is 
the infinity of statistical moments. Indeed, the identity oiU + U' and U + U in the sense 
of moments is possible only if the variance of U is infinite. This means that probability 
density has a long power-law tail. 

In practice one hopes that random rates are somehow bounded and have finite cu- 
mulants and moments. From the other hand, even if a statistical ensemble has infinite 
moments, this infinity never could be observed, at least while finite time intervals. Any 
moment practically determined with the help of time averaging always is finite with unit 
probability. 

But no theory can do without ensemble language. Therefore, we need in an ensemble 
which would be close to the Cauchy ensemble but possess finite moments and be applicable 
to real time-limited observations. 

Clearly, such an ensemble will allow for as large AQ(t) deviations from its most prob- 
able behaviour as long is the observation time, in accordance with the temporal scaling, 
AQ(t) oc t , of the ideal Cauchy random walk. Hence, this ensemble inevitably includes 
long-living statistical correlations although thats express nothing mote than the perfect 
randomness of the process. 

E. SELF-FORGETTING AND SCALE INVARIANCE 

The universal features of such statistics can be deduced from the simple requirements 
to follow [6,8]. 

a) . The ensemble average value of the accumulated amount is strictly proportional to 
time, (AQ(t)) = U t , that is the mean value of time-averaged rate does not depend on t 
, (U(t)) = U . 

b) . In the limit t/ro — > oo , the quantities AQ and U are constituted by many causally 
independent incidents, so the probability distribution of AQ should acquire an infinitely 
divisible form. In accordance with the general Levy-Khinchin representation of infinitely 
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divisible distributions (for one-directed asymmetrical walks with AQ > ) [12], the CF 
of AQ can be written as 



The essence of this formula is that function G is non-negative, G(q, t) > . Besides, 
because of condition a) , the equality / °° G(q,t)dq = 1 should be satisfied. 

c) . Due to t » t , AQ(t) is believed to asymptotically behave as a random walk with 
independent increments. This means that logarithm of CF is approximately proportional 
to time, In (exp{ikAQ}) oc t , i.e. G(q,t) becomes almost time independent at large t . 

d) . AQ(t) should asymptotically behave in a scale- invariant way, because of indiffer- 
ence to the short-time scale To and absence of larger characteristic scales. The concrete 
form of scale invariance is prompted by the mean law (AQ) oc t and looks as 



where symbol ~ denotes the approximate statistical equivalence. Thus though being 
random AQ(t) grows nearly proportionally to time. 

The physical meaning of this assertion was already discussed: a process what is cur- 
rently forgetting its previous accumulations (number of events, Brownian displacement, 
etc.) is unable to make difference between its mean part and deviation part. Any result 
of previous behavior is equally acceptable, regardless of its difference from ensemble expec- 
tations! Therefore it is more correct to say that just the average (AQ(t)) oc t reflects the 
scaling AQ(t) oc t of random increments. That is why we must not especially carry out 
ik (AQ) from In (exp{ikAQ}) . 

e). In the limit t/r — > oo the time-averaged rate tends to U . Thus we even want 
the process be ergodic. We will see that the ergodicity automatically follows from a)-d) , 
but it is satisfied only in a logarithmically weak sense, as well as c) and d) do. 



If being formally punctual the scale-invariance condition d) would mean that 
XG(q,\t) = G(q/\,t), that is G(q, t) = q~ 1 ^(q/Vt) . Here $ is formally arbitrary 
function and V is some characteristical rate introduced to make the argument q/Vt di- 
mensionless. But this purely invariant form is incompatible with conditions b) and c). 
Indeed, if $(0) ^ then the integration diverges at q — > , while if $(0) = then the 
condition c) will not be satisfied. 

Hence we should allow for a slight violation of the scale invariance. We can expect that 
necessary violation will be really weak because the divergency is logarithmically weak. 

The correction must cut the divergency at small values of q , q < Qq , with Q being 
some characteristical "microscopic" scale for measuring AQ . Such formal correction 
reflects a really inevitable violation of the invariance at small increments < Qq as well 
as at small time intervals. For instance, one could simply replace the integration from 
zero to infinity by integration starting at q = Q (there are many other possibilities but 
results. will be practically the same). Then condition a) yields 




AQ(Xt) ~ XAQ(t) 



F. CHARACTERISTIC FUNCTION OF 
PURELY RANDOM FLOW OF EVENTS 
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1 = r -*(£i)dq « HQo/Vt) ln(^) » $(0) ln(-) 

where r = <5o/^ , an d the inequality i >> r is supposed. 

From here it follows that the correction is forced to include additional logarithmical 
dependence of $(0) and so of G(q,t) on time. This dependence can be separated as the 
multiplier 

A(t) « [hi-]- 1 

TO 

if make the change =>■ A(t)^(y^) accompanied by the requirement $(0) = 1 . 

If all the cumulants and statistical moments of AQ are finite then should decrease 
from unit to zero in a sufficiently fast way. Details of this function are of no principal 
importance. For example, one can choose = exp(— x) . Finally, the simple estimates 
of the integral give 

(exp{ikAQ}) m exp{-ikU tA(t) ln(— - ikQ )} = 

is 

= exp{ikU t - ikU tA(t) ln(l - ikVt)} 
under the conditions t » er and \k\Q « 1 , with the latter corresponding to AQ » 

Qo- 

Evidently, this is the CF of probability distribution which is as much close to ideal 
Cauchy distribution as the limitations of statistical moments do allow. There is also 
slightly other way to avoid the divergency [8], which results in the CF {exp{ikAQ}) m 
exp{ikU t[l — ikVt]~ A( -^} . The corresponding probability distribution resembles Levy- 
Khinchin stable distributions with Levy's alpha-parameter close to 1, = 1 — A(t) , 
and practically reduces to the above Cauchy type distribution if A(t) « 1 . One more 
possible modification was obtained in [10]. 

G. LONG-RANGE STATISTICS OF PURELY 
RANDOM FLOW OF EVENTS 

The time scale r should be treated just as life-time of the memory of random walk: at 
t » t the transition to almost scale invariant behaviour takes place. At once, the factor 
^ in logarithm serves as the cut parameter. If we neglected it the CF would transform 
into 

(exp{ikAQ}) ~ exp{-ikU tA(t) \n(-ikQ )} 

what corresponds to pure Cauchy statistics with infinite moments. 

The probabilistic properties of AQ significantly depend on the parameter K,(t) = 
V/[U A(t)] = Q /[T U A(t)} . If k « 1 then the distribution is almost Gaussian. In 
opposite case, k » 1 , it approaches to ideal asymmetric Cauchy distribution and its 
density has the long tail U tA(t)AQ~ 2 although cut at far AQ values [6,8]. In this 
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case all the probabilistic properties become almost insensitive to how far the long tail 
is cut. For example, the estimates take place P{U > Uq} ~ - arctan(7r/ ln/t) , and 
P{U > 2U Q } A(t) exp(—U T /Q ) where P{} denotes the probability. 

The essential feature of the statistics is a difference between the mean and most probable 
values of AQ as well of the rate U = AQ/t . In particular, at k > 1 this difference is 
decsribed by U — U mp ps U A(t) In n(t) , with U mp being the most probable rate. 

H. CAUCHY STATISTICS AND 1/F-NOISE 

In contrary to probabilities, the second and higher statistical moments are exteremely 
sensitive to the long tail being determined mainly by bad untypical occurences. From the 
obtained CF one gets 

(AQ 2 ) - (AQ) 2 ^ 2U % 2 A(t) , (U 2 ) - (U) 2 ~ 2C/ — A(t) (t » r„) 

Therefore, the variance of time-averaged rate U (t) tends to zero but only in logarith- 
mically slow way. This means that the rate undergoes low-frequency 1/f fluctuations. It 
is easy to show that their relative spectrum is Ssu(f) ~ (<3o/^o r o)/ _1 ( m 2nfT )~ 2 , at 
fr « 1 • 

The only free dimensionless parameter of the statistics is Qo/Uqt . It determines 
both the 1/f-noise level and the degree of non-Gaussianity of the statistics. Therefore, 
the close relations between probabilistic properties and spectral properties of memoryless 
random flow of events do exist, and 1/f-noise is as strong as non-Gaussian. Both the 
1/f spectrum and bad statistics mutually come from the absence of long-living memory. 
Their quantitative connection is the manifestation of highly irregular fractal structure of 
realizations of AQ(t) [8] similar to that of Levy flights with a L ps 1 . 

I. ABSOLUTE TIME SCALES AND SCALELESS FLUCTUATIONS 

In general, any random walk subjected to a stable statistics [12] possesses self-similar 
fractal structure. But, unfortunately, a pure mathematical Levy flight has no intrinsic 
time scale at all, and so its rate looks as a white noise. Besides, at alpha-parameter 
a L < 1 even its mean value turns into infinity. This is too idealized model to apply it 
to real processes which always have some absolute lower time scale. We showed that 
in the case ai = 1 — the existence of such a scale automatically leads to appearance 
of logarithmically time-dependent factors like A(t) responsible for long-living scaleless 
statistical correlations. One comes to the idealized statistics only in the limit r — > , 
U oc 1/A(t) -> oo . 

J. LOW-FREQUENCY NOISE AS DETERMINED 
BY SHORT-TIME SCALES 

Formally Qo/UqTq is arbitrary parameter, but in concrete applications it is possible 
to suggest some natural estimates for it. For example, let AQ(t) be formed by well 
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distnguishable identical discrete events (dozes, electron jumps, etc.). In this situation it 
is natural to choose Qq = 1 , and Uq ~ l/r c , with r c being the mean time separation of 
events, that is the correlation time of the Poissonian short-range statistics of u(t) . Then 
Ssu(f) ~ ( T c/To)/ _1 (hi27r/r )~ 2 , and the relative 1/f- noise level is determined by ratio 
of mean distance between events and the memory duration, i.e. by the amount, ~ t /t c 
, of recent events which could be kept in memory. 

The corresponding mean square relative uncertainty of the rate if measured during 
t » t is (5U 2 ) ~ (j c /tq) (In t/r ) _1 . Hence, the rate can not be predicted with much 

better precision ( « \J r c j r ) than what is available for the system itself and achievable 
by its self- observation during its memory life-time. 

The particular case t t c corresponds to minimally short memory, with low- 
frequency noise Ssu(f) ~ j > where a = (In27r/Y )~ 2 , a ~ 0.02 at fr ~ 10 -4 and 
a ~ 0.002 at fr ~ 10~ 9 -j- 10~ 13 (what is typical for electric noises). It can be shown 
that at a > 0.0001 the probability density of U has noticable long tail. 

The constructed model is able to give satisfactory quantitative estimates [6-8] of both 
the spectrum of 1/f-noise and its statistics which can be experimentally investigated [22- 
24] in small carbon resistors, thin metal whiskers and some other systems (in particular, 
a shape of experimental probabilistic histograms and the difference between average and 
most probable values can be correctly described). 

K. 1/F-NOISE IN DNA 

Now we are able to hypothetical^/ interprete also the example ix) [17,18]. The 1/f- 
noise in a distribution of any given sort of bases over DNA sequence could mean merely 
that only more or less localized pieces of DNA serve for encoding an information, while 
the amounts of bases accumulated by long pieces and total number of given bases in DNA 
are of no importance. 

If this is true, we concern only statistical long correlations which themselves do not 
have any specific significance. But, in accordance with above phenomenology, the mag- 
nitude of 1/f-noise could contain an information about typical length of actually self- 
connected pieces. 

In any case, the difference between statistical and actual (causal, structural) correla- 
tions should be taken into account. The paper [18] is one a very few works we know where 
this difference is realized. 

L. COIN TOSSING AND LASER BEAM 

Of course, the statistics under consideration is only the principal possibilty, not a 
necessity. Everybody knows about random flows of events which keep no place for 1/f- 
noise. For instance, the mathematical coin tossing first investigated by J.Bernoulli in his 
"Ars Conjectandi" in 1713. 

However, one could say nothing about flipping a physical coin tossing, if do not assume 
that the events are "statisticaly independent" and the probabilities of heads are one half, 
thus coming back to the mathematical model. May be, the probabilities differ from one 
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half, but this is of no importance if they are thought as purely certain. The essence of the 
"statistical independency" is that one identifies the probabilities of heads with the "the 
probabilities per unit time" (i.e. per unit tossing) which eventually coinside with limits 
of averaging over many tossings. 

Therefore, to say that events are statistically independent is the same as to assume that 
these limits exist and besides are achieved in sufficiently fast way (not in logarithmical 
one). Hence, the only factual Bernoulli's achievement was the law what describes the 
approach to the limit, provided it does exist, not the proof of its existence. 

What is for the physical coin, may be, in this case the flicker floor is really negligible. 
But the opposite reality does not imply that there are some slow processes, it means 
rather that the ideal certainty of the probabilities materialized in ideal constancy of the 
coin takes no place. 

In fact, the ideal constancy of coin serves as infinitely long-living causal correlation 
what just ensures the absence of 1/f-noise in the model. The Poissonian statistics of 
photons in an ideal laser beam is the more modern example of Bernoullian flow of events. 
The absence of 1/f-noise in photon flow in ideal beam also is due to hidden infinitely 
long causal correlation, namely, to what all photons belong to one and the same quantum 
state. Similarly, an electric supercurrent has no 1/f fluctuations, in contrary to normal 
current in superconductors. 

M. QUANTUM DECAY 

From the viewpoint of conventional statistical interpretation of quantum mechanics, 
the random flow of decays mentioned in the example a) is very similar to coin toss- 
ing. One needs in definite probability of decay per unit time and, besides, in statistical 
independency of decays in order to predict the Poissonian statistics of decays. The 1/f 
fluctuations (observed at frequencies below 10 _3 -=-10 -4 Hz ) are in principal contradiction 
to this theoretical scheme. So we can conclude that the hypothesis about the statistical 
independency is wrong. In other words, the decay probability has no absolutely definite 
value. The observed 1/f-noise is quite equivalent to 1/f fluctuations of probability of 
decay. 

In fact, one lacks the statistical independency in order to connect the decay probability 
calculated in quantum mechanics and imaginary ensemble of identical experiments. But 
if we wanted to prove the statistical independency of two events we would need in the 
more complicated assumption that any two pairs of events are statistically independent, 
and so on up to infinity. Therefore the statistical independency can not be grounded in a 
logical way. This circumstance is noted in the book [25] . 

It is the open question whether an uncompleteness of quantum mechanics or primitiv- 
ity of decay models is responsible for the discrepancy. Another problem was pointed out 
by Prigogine [26] . It is well known that time dependence of the decay probability includes 
non-exponential (power-law) tail. It relates to times much larger than the life-time and 
so it itself is practically unobservable. Nevertheless, it is inconsistent with the princi- 
ple that quantum decay does not remember the past: if one knows that a nucleus has 
not yet decayed then any expectations for future should not depend on how long it was 
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being safe. As a consequnce, the probability of no decay should be purely exponential, 
P(t) = exp(-Ti) . 

In our opinion, it may be possible to avoid at once both the contradictions in the 
unified way if suppose (or better deduce from quantum mechanics) that the probability 
of decay per unit time, 7 , is not purely certain, and that its uncertainty is just as large 
as necessary in order to disguise the nonexponential tail. In other words, we should make 
difference between the most probable behaviour of experimentally registered flow of decays 
and its ensemble-averaged behaviour, as in general in memoryless flows of events. Then 
the quantitative measure of this difference is nothing but that of 1/f-noise. 

N. EQUILIBRIUM THERMAL 1/F-NOISE 

The equilibrium voltage noise in a disconnected tunnel junction is constituted by 
random jumps of electrons from left to right-hand side of the junction and in opposite 
direction. Let N + (t) and A_(t) be the numbers of these jumps during some time interval 
t . If N + > AL then the electric voltage arises between sides, u = e(N + — NJ)/C , with C 
being their electric capacity, which makes one sort of jumps more probable than another. 
Due to this thermodynamical opposite reaction, after a characteristic relaxation time of 
order of r c = RC , with R being the junction resistance, the system returns to a vicinity 
of the chargeless state. Therefore the difference N + — AL always is kept near its average 
zero value, and the time-smoothed current e(A+ — NJ)/t turns to zero at t » r c . 

But let it happened that the sum N(t) = N + + AL ps 2AL essentially deviates from 
its average value, while the difference remains close to zero. Such a fluctuation causes 
neither a charging nor any other thermodynamical consequences. All stay as before! 
So we can expect that later this fluctuation will not be "compensated and damped" by 
some oppositely aimed process. Therefore, in reality the time-smoothed value N(t)/t has 
no reasons to tend to its "theoretical" ensemble averaged value, moreover, it have no 
explicitly predictable limit at all [10]. 

Now notice that both types of jumps give an equal positive contribution to the noise 
intensity. The latter is determined by the summary number of jumps per unit time, 
S = e 2 N/t = e 2 (AL + AL)/t , with S being the instantaneous (time-averaged) power 
spectrum. Hence, we come to white noise but whose intensity undergoes low-frequency 
1/f-fluctuations, i.e. to flicker "noise of noise". 

O. 1/F-NOISE AS UNCERTAINTY 
OF KINETICAL VALUES 

It is important to underline that in all the above examples the values subjected to 
1 / f-fluctuations could be qualified as kinetical values. This is principal peculiarity of 1/f- 
noise without slow processes. In contrary to dynamical and thermodynamical variables, 
kinetical values characterize the statistics of random transitions between various states of 
the system, not statistics of the states themselves. 

In general, good ergodicity with respect to visiting the states does not imply the same 
with respect to transitions. Even in the simplest case, when a model distinquishes only 
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two states and each of them occupies well certain portion of time, the rate of random 
jumps between the states may be uncertain. 

None kinetical value could be referred to a definite time moment. The matter is that 
any kinetical value takes a shape only during a finite time interval longer than memory 
life-times. This circumstance means that kinetical values have no their own characteristic 
time scales. As a consequence, their fluctuations naturally acquire a scale invariance with 
a scaleless power-law low-frequency spectra. 

P. 1/F-NOISE IN NON-EQUILIBRIUM STATE 

In thermodynamical systems every sort of kinetical events (random transitions from 
one state to another) is accompanied by time-reversed events, and a freedom of random 
behavior is restricted by definite ratio of amounts of opposite events (by detailed balance, 
in equilibrium case). This restriction arises from reversibility of underlying dynamics, 
always are under thermodynamical control and can be taken into account with the help 
of fluctuation-dissipation relations. The example is the ratio N-/N+ ~ 1 in equilibrium 
electric junction. 

But fluctuations in the summary rate of direct and reversed kinetical events does not 
destroy their balance and so does not cause complete counteraction. 

If an external voltage U is applied to the junction, then undumped fluctuations in the 
total rate of jumps result in not only noise of noise, but also strictly in the voltage and 
current. Now the relation AL ~ N + exp(— eU/T) is thermodynamically kept, instead of 
N_ ~ N + , and the current, e(N + - N_)/t = etfmh(eU/2T)N/t = tanh(eU/2T)S/e oc 
S , begins to reflect uncertainty of the total number of jumps per unit time. In this 
nonequilibrium state the voltage spectrum acquires 1/f contribution in addition to the 
white background, and 1/f-noise can be viewed directly in the two-point correlators as 
effective resistance fluctuations. 

IV. BROWNIAN MOTION REVISED 

A. QUAZI-GAUSSIAN RANDOM WALKS 

Let u(t) be equilibrium noise short-correlated in the sense of two-point correlators. 
Then the integral Q(t) = $lu(t)dt is symmetrical random walk like Brownian motion. Its 
characteristic fractal properties are determined by the scaling Q(t) oc \ft , at least if t 
sufficiently exceeds tq, with tq being maximal correlation time of u(t). As it was shown 
in [4,5], if only one excludes very formal case of ideally Gaussian u(t) , then this scaling 
results in 1/f-noise of instantaneous power spectrum, S(t) , of the noise u(t) and in definite 
non-Gaussian long-range statistics of Q(t) (termed quazi- Gaussian) which resembles the 
statistics of Levy flights with alpha-paremeter slightly smaller than 2, a L = 2 — . One of 
many examples is the equilibrium voltage noise in disconnected electric junction. Another 
example is the equilibrium Brownian motion, with u(t), Q(t) and D(t) = S(t)/2 being 
the velocity, displacement of a Brownian particle and its current diffusivity, respectively. 
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Principally, the statistics of Q(t), like that of the above analysed one-directed random 
walk, is derivable from rather trivial requirements to be satisfied asymptotically, at t » 
To- 

a) The standard diffusion law takes place, (Q 2 (t)) = 2D t . 

b) Q(t) behaves as infinitely divisible random value. 

c) It approximately behaves as random walk with independent increments. 

d) Due to forgetting the past the system does not distinguish what is a proper rate of 
the walk and what is deviation from it. Any currently happened location of the walk is 
equally permissible as starting point for further adventures! Therefore, the scaling Q(t) oc 
\/t is peculiar for the statistics as a whole (not only for the mean square displacement). 
This can be expressed in the form of approximate statistical identity 

Q(Xt) ~ V\Q(t) 

These requirements imply nearly stable probability distribution whose characteristic 
function (CF) looks as [5,7,8] 

(exp{ikQ}} « exp{D k 2 tA(t) ln( — + k 2 q 2 )} = 

= exp{-D k 2 t + D k 2 tA(t) ln(l + D'k 2 t)} 

Here A(t) ~ l/ln(t/r ) , t » r , k 2 q 2 « 1 , D' — q 2 /r , and the scales q and r 
serve as lower spatial and temporal bounds of the scale invariance of the random walk. 
Under the formal limit r — > , D oc \/A(t) — > oo , this CF turns into CF of pure quazi- 
Gaussian distribution, exp{— Dtk 2 \w.{l/k 2 ql)} , with infinite second and higher even-order 
moments and without characteristic time at all. This random walk constructed in [4,5] is 
absent in standard classification of Levy flights [12]. Independently on us it was found in 
the mathematical work [27]. But many of applications need in random walk models which 
possess finite statistical moments and non-zero absolute lower time scale with consequently 
deformed scale invariance, as in the model under consideration. 

B. FOUR-POINT CUMULANTS 
AND DIFFUSIVITY FLUCTUATIONS 

The main information about 1/f- noise is contained in unusually behaving 4-order cu- 
mulant of the displacement what follows from the above CF, 

(Q {A \t)) = (Q(t),Q(t),Q(t),Q(t)) Kt 2 A(t) « tXt/To)]" 1 

In general, the standard variance (Q 2 (t)) = ID^t can be accompanied by non-standard 
4-point cumulant like (^Q^{t)^ oc t l+1 with 7 > 0. At any 7 > such the asymptotics is 
formally equivalent to low-frequency fluctuations of instantaneous power spectrum (dif- 
fusivity) of u(t) . Indeed, by the formal definition of cumulants, 

(exp[ikQ(t)]) = exp[-k 2 (Q 2 (t)) /2 + k A (Q (4) (t)) /24 - ...] 
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From the other hand, one can phenomenologically treat u(t) as a shot noise what would 
be Gaussian if its power spectrum was constant, S(t) = 2D , but becomes non-Gaussian 
because of its relatively slow random modulations, S(t) = 2D(t). Correspondingly, the 
CF can be represented by double averaging, first over fast Gaussian noise under fixed 
realization D(t) and then over random D(t) : 

(exp[ifcQ(t)]> ~ ((exp[-k 2 J* D{t')dt']j^ = 



= exp[-k 2 D t + k 4 f f K D (ti - t 2 )dt 1 dt 2 /2 - ...] 
Jo Jo 

Here the doubled brackets denote the second stage of averaging, D = (D(t)) , the 
inequality t » r c is implied, and K D {r) is the effective correlation function of fluctuating 
diffusivity D(t) , 

Ko{t) = ^( g(4)(t) ) = 111 11 «t)Mr')Mr"),u(0))dr'dr" 

The unusual 4-point cumulant means slow asymptotics K D {r) oc r 7_1 and thus /~ 7 
low-frequency spectrum of diffusivity fluctuations. The case 7 > 1 formally corresponds to 
non-integrable spectrum like that of nonstationary random process, but this is only imag- 
inary non-stationarity because u(t) is purely stationary random process (all its cumulants 
possesses time translation invariance). 



C. MICROSCOPIC TIME SCALES 
AND MAGNITUDE OF 1/F-NOISE 

In case of quazi- Gaussian random walk 7 — > 1 , K d {t) oc A(t) , and the low- frequency 
part of spectrum of relative diffusivity fluctuations is S$£>(f) ~ a(f)/f , where a(f) = 
(g 7A>To)[ln27r/T ]- 2 . 

Let us estimate this noise taking in mind literally Brownian motion of a microscopic 
particle. The mean diffusivity can be expressed as D ss q 2 /r c with q c and r c being its 
mean free path and free flight time, respectively. Clearly, r c coinsides with correlation 
time of the velocity u(t). It is quite reasonable to suppose that spatial scale invariance 
takes beginning just from one free path, that is q ~ q c . Then magnitude of 1/f- noise is 
determined by the ratio r c /r . 

Again we see that 1/f- noise is as strong as long is the memory duration as compared 
characteristic time separation of elementary kinetical events (interactions, collisions, etc.). 
The smaller is a number of previous kinetical events remembered and taken under control 
by a mechanism what produces white noise the greater is the accompanying 1/f- noise. 
In fact, one can be convinced of validity of this statement [7] with respect to different 
particular cases of Brownian motion of charge carriers in semiconductors and metals, to 
charge transport in normal junctions and superconducting Josephson junctions, and to 
Brownian particle in liquids (in the latter cases the difference between q and q c can be 
significant too). 
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The special case, tq ps t c , when scale invariance starts immediately after one free flight 
time, corresponds to natural basic level of diffusivity noise. This is the case for electrons 
in a slightly doped semiconductor, due to inelasticity of collisions. Taking into account 
that r c ~ 1CT 12 s , at / ~ 1 Hz (as well as at frequencies a few orders larger or smaller) 
one finds SsD(f) ~ 0.0015// . This estimate [4,5] is in good agreement with empirical 
Hooge constant a m 0.002 [1,2]. 

D. PROBABILITY DISTRIBUTION 

OF BROWNIAN MOTION 

In contrary to 4-th and higher cumulants, the probability density W(Q,t) of quazi- 
Gaussian random walk always has only small quantitative differences from that of ideal 
Brownian motion, W G (Q,t) = (ArrDt)- 1 / 2 exp(-Q 2 /4Dt). If D' > D then W(Q,t) has 
noticable power-law tails whose extent approximately equals to \/2D't but magnitude 
does not depend on D'/Dq , i.e. on the 1/f-noise level. In the limit D'/Dq — > oo, for 
the region Q 2 » (Q 2 ), one can get W(Q,t) « W G (Q,t) + 2D tA(t)/ \Q\ 3 . From here 
the probability P of that Q 2 is three times larger than (Q 2 ) = 2D t can be estimated as 
P < Pq + [3 ln(t/r )] _1 , where P G means similar probability for ideal Brownian walk, 
P G ~ 0.1 . Therefore, at t/r > 10 3 the difference between P and P G never could 
exceed 0.05. Analogouosly, one can find that at t/r > 10 3 always P{\Q\ > 4:^/2Doi} w 
[Wlnt/ro]- 1 < 0.01 . 

E. STATISTICAL SENSE OF LONG 
CORRELATIONS OF KINETICAL VALUES 

Hence, the probability of noticably far deviations from standard behaviour does not 
exceed one percent even in the worst case (infinite 1/f-noise). Nevertheless, just such 
rather unprobable occasions determine the asymptotics of Kc(t) and the level of diffusivity 
1/f-noise. 

Evidenly, the long tail K D (t) oc A(t) reflects only the statistical weight of untypical 
realizations of the random walk, not some actual correlations inside realizations. More 
explicitly, any concrete realization looks as sometimes typical and sometimes rare, and 
because of its fractal structure every variant of its behaviour during longer time interval 
is constituted by both typical and rare contributions from shorter subintervals. 

The same can be said also about fluctuations of other kinetical values. From the 
viewpoint of underlying dynamics, their correlation functions are roughly phenomenolog- 
ical characteristics. However, one can present rigorous definitions for these correlators in 
terms of four-point (or higher order) cumulants (or related ensemble averages) of purely 
dynamical variables. 

V. DYNAMICS AND 1/F-NOISE 
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A. SOME QUESTIONS 



Of course, we are highly interested in 1/f- noise derived from Hamiltonian dynamics. 
Unfortunately, none of realistic Hamiltonian models allows for explicite calculations of 
4-point cumulants even for Hibbsian equilibrium ensemble. Moreover, this problem never 
was under careful consideration. Instead, everybody wanted to somehow reduce the statis- 
tics to two-point correlators thus diminishing specific contribution of 4-point cumulants. 
Various means were attracted to this aim, either ansatzes like molecular chaos, random 
phases, thermodynamical limit, continuous spectra of energy levels, diagonal singularity 
and others, or exclusion of some interactions what do not influence low-order correlations 
under interest. 

Some of well investigated low-dimensional nonlinear mappings have the intermittency 
regimes accompanied by 1/f- noise, but this noise is connected with long periods of regular 
motion, that is with long memory and slow processes (for example, see [28]). Therefore 
it is not surprising that this 1/f-noise becomes destroyed (saturated at zero frequency) 
under influence of external white noise. So it is not the 1/f-noise we are interested in 
whose low-frequency behaviour should not be supressed by noisy perturbation. 

At present the axiomatic theory of the so-called Anosov systems is developed [29,30] 
based on works by Sinay on scattering billiards and ideas by Anosov. These systems are 
characterized by exponential instability and divergency of phase trajectories and their 
mixing in phase space, possibly in company with time reversibility and phase volume con- 
servation (what is most interesting for us). In such systems 1/f-noise could be connected 
with random temporal non-uniformity of the rate of exponential instability. As far as we 
know, this hypothetical possibility was not under investigation. In contrary, the attention 
was concentrated on the systems and phase space divisions which produce Bernoullian 
flows of events. 

However, as it is was shown by Ornstein [13], in principle generalized random flows 
of events (K-flows) do exist which can not be recoded into Bernoullian flows. One could 
search for 1/f-noise in concrete realizations of such generalized flows, but this subject is 
still apart of achievements of the theory. 

B. WHY KINETICS LOST 1/F-NOISE ? 

Because none realistic dynamics can be exactly solved, physicists resort to kinetical 
models of relaxation, noise and irreversible phenomena. Everybody knows how this is 
carried out. First one chooses a reasonable coarsened characterization of the system, 
secondly, analyses what variants of its evolution may occur during some short time step At 
if start from an uncompletely described state, then evaluates probabilities of the variants 
and finally keeps these probabilities as governing parameters for a long-time evolution 
during many steps. 

This theoretical scheme was expressively criticized by M.Kac [31]. Clearly, its main 
defect is its last ansatz, because it means that probability of any chain of transitions 
between states is equalized to product of probabilities of marginal transitions, i.e. that 
successive transitions are considered as statistically independent events. One could say 
that such a model divides the real trajectories to short pieces and randomly permutes 
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these pieces at each time step At. As a consequence, the model can correctly describe 
short correlations and relaxation, but in general is unable to correctly reproduce a long 
time behavior of trajectories. 

Usually one sees nothing criminal in this simplification. But we insist that long-living 
statistical dependencies and non-Bernoullian behaviour may mean merely that trajec- 
tories take liberties just because of no really long-living memory. Then the statistical 
independency looks as too rough approximation. Eventually, this approximation results 
in the loss of 1/f- noise. 

We see that modern theory is not kind with respect to 1/f-noise. Nevertheless, we 
have pointed out two Hamiltonian systems, both being basic models of statistical me- 
chanics, which produce 1/f-fluctuations of kinetical values, namely, a slightly non-ideal 
gas [9] and slightly unharmonic phononic system [11]. In our opinion, every many-particle 
Hamiltonian model which really bears a relaxation inevitably bears also 1/f fluctuations 
of the rate of relaxation. 



Many years ago M.Kac [31] suggested the model which represents a good toy analogy 
of the problems always accompanying derivation of relaxation, noise and kinetics from 
dynamics. The system is the ring chain of n » 1 boxes each containing the only black 
or white ball. At every moment of discrete time t — .. — 1,0, 1,2. ..the balls displace to 
left along the ring. The boxes are characterized by time- independent parameters Oj = 
±1 . If Oj — — 1 then j-th box inverts the colour of any ball what leaves it. The 
dynamical variables £j(t) = ±1 are defined so that £j(t) = 1 if j-th box contains a white 
ball at time moment t and = —1 if it contains a black ball. The equations of 

motion are easily soluble: £j(t + 1) = o\,_i£j_i(i) (the number j = —1 is equivalent to 
j = n — 1). The equilibrium statistical ensemble is defined by conditions that £j are 
statistically independent if considered at one and the same time moment and £j = ±1 
with equal probabilities 1/2 , and that constant parameters Oj equals to 1 or —1 with 
probabilities p, 1 — p also being independent one on another. Evidently, the dynamics is 
time-reversible and the ensemble is time- invariant. We confine ourselves by the simplest 
case p — 1/2 . 

The quantity under macroscopic observation is the mixed colour u(t) = [J2j £j{t)]/ V™ 
. It is easy to verify that (u(t)) = and that from the point of view of second-order 
correlators u(t) represents a white noise: (u(ti)u(t2)) = S tl t 2 , with Kronecker symbol on 
the right-hand side. 

Kac demonstrated [31] that this model allows for wide analogies with serious statistical 
mechanics and for illustrations of approximate methods usually attracted when building 
kinetics. But Kac did not note that u(t) is very far from white noise in the sense of 
higher-order statistics. 

Let us consider the discrete random walk Q(t) = Y?k=t +i u (k) at 1 < t < n. Obviously, 
(Q 2 (t)) = 2D t, with diffusivity D = 1/2. If this process was asymptotically gaussian 
and subjected to the central limit theorem at t » 1, the four-order cumulant of Q(t) 
would grow also as a linear function of time, (Q^(t)) oc t. For instance, if u(t) was 
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purely white noise, one would get (^Q^(t)J = —2t/n . However, the explicit calculation 
lead to 



(Q^(t)) = (Q 4 (t)) - 3 (Q 2 (t)) 2 = 2{2t 2 - l)(t - 2)/n oc t 3 (t » 1) 

(this easy reproducable calculation was performed ten years ago and sent to 1988 Confer- 
ence on noise in physical systems in Vilnius, but, unfortunately, was not accepted either 
there or for journal publications). 

This result proves the existence of infinitely long-living four-point correlations (cumu- 
lants) of the noise u(t), and means that standard kinetical approaches can not present 
a qualitatively correct approximation (let being quantitatively good) of the noise. From 
the viewpoint of external observer, such 4-th cumulant looks as manifestation of low- 
frequency fluctuations of diffusivity, i.e. of instantaneous power spectrum of u(t). In 
accordance with the above formulas, the corresponding spectrum is oc f~ 2 . The factor 
\jn is analogous to what takes place in the spectrum of relative resistance 1/f-noise in 
bulk conductors (example i) ). 

D. NON-STATIONARY NOISE 
IN STATIONARY SYSTEM 

Perhaply, this noise could not be qualified as what is due to forgetting the past, sooner 
the actual correlations carried in by random parameters <Jj do work. Nevertheless, the 
obtained l// 2 -noise of noise is principally important by two reasons. First, it gives the 
example of completely calculatable long-living higher-order correlations in time-reversible 
dynamics. Secondly, it helps to understand why seemingly non-stationary fluctuations of 
kinetical values can occur even in equilibrium system. 

With no doubts, the non-integrable spectrum 1/f 2 belongs to stationary system . 
Indeed, the statistics of Q(t) does not depend on the arbitrary initial time moment to- 
The spectrum 1/f 2 reflects only the role of duration t in the measurement of kinetical 
value (diffusivity) and means that the longer is the measurement the worse is its result. 
The similar strange property characterizes Levy flights with q £ < 1 [12]. 

Of course, the guessing is that any kinetical value is delocalised in time. Therefore 
its bad behaviour does not contradict to ergodicity of visits of states. Moreover, we 
can expect that at fixed duration t the four-point cumulant calculated with the help of 
ensemble averaging coinsides with that obtained from the time smoothing over t . 

Besides, one should remember that non-stationarity in the sense of statistical moments 
coexists with stationarity of the most probable behavior. But because of non-stationarity 
the result of experiments can be sensitive to details of the measuring procedure [8,9]. 

E. 1/F-NOISE OF FRICTION AND OF LIGHT 
SCATTERING IN PHONONIC SYSTEM 

A quartz crystal is the system where thermodynamically equlibrium and nonequilib- 
rium manifestations of one and the same 1/f-noise are observed by different methods. 
As it is known for a long time [14,15], when quartz is used as frequency stabilizator in 
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electronic devices, its quality is restricted by 1/ f-fluctuations of inner mechanical friction 
and consequently of dissipation. Not long ago the 1 /f-fluctuations of intensity of the laser 
light scattered by equilibrium quartz were discovered [16]. 

In [11] we showed that the statistical correlators which describe the fluctuations of 
dissipation can be generally connected with four-point equilibrium cumulants which de- 
scribe thermal fluctuations of quantum probabilities of photon scattering by phonons. It 
was shown also that both the friction and light scattering 1 /f-fluctuations can be reduced 
to thats of the relaxation rates (life-times) of phononic modes. In their turn, this fluctu- 
ations originate from the same phonon interactions what produce exponential instability, 
continuous spectra and the relaxation phenomenon itself. 

The logical way to 1/f-noise looks as follows. Due to exponential instability, the in- 
teraction of given test phononic mode with the thermal bath consisting of other modes 
produces both its noisy disturbation and damping. In zero approximation, the friction is 
referred to linear response of the bath to test mode and is non-random, and correspond- 
ingly the noise is Gaussian. The linear response seems good approximation because the 
influence of the only test mode is l/y/N times smaller than summary interactions inside 
the bath, with N being total number of modes, N ~ 10 20 . 

Next, the equations which govern the linear response look as equations of coupled har- 
monic oscillators without friction but with fluctuating couplings. Even being Gaussian 
this fluctuations cause a stochastic instability of the response, in the sense of second and 
higher statistical momemts. Therefore in better approximation the friction becomes ran- 
dom with infinite variance and besides with infinitely long non-locality. As a consequence, 
the infinitely long-living four-point statistical correlations of phononic variables do occur. 

In a more developed picture, one should consider the completely nonlinear response. 
But, as it was argued in [11], the long-living many-point correlations survive and corre- 
spond to giant low-frequency fluctuations of phase diffusivities of phonon modes and thus 
of their relaxation rates (relaxation times). 

VI. GENERALIZED FLUCTUATION- 
-DISSIPATION RELATIONS 

A. PROBABILISTIC RELATIONS 

In any theory concerning 1/f-noise in thermodynamical systems the exact generalized 
fluctuation-dissipation relations (FDR) can be useful which were considered in [32-34] 
(see also [7,9,11,38]). Except classical fluctuation-dissipation theorem, Kubo formulas and 
Onsager relations, FDR include infinitely many additional connections between nonlinear 
responces and higher-order equilibrium and nonequilibrium cumulants. 

All the variety of FDR can be expressed with only relation what follows from the 
phase volume conservation and time reversibility of Hamiltonian dynamics. For instance, 
in case of initially equilibrium Hibbsian canonical ensemble later perturbed by some time- 
dependent variation of the Hamiltonian, the producing FDR can be written [32-34] as 

U{Trajectory; Forces} exp(—Work/T) = U{Trajectory; Forces} 
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There U{Trajectory;Forces} is the probability of realization of some process 
(Trajectory) under given external Forces=Perturbations+Conditions (in the sense of 
probability density in functional space), Work is the work produced by the perturbation 
during this process depending on Trajectory and Forces , T is the temperature of initial 
Hibbsian probability distribution, and Yl{Trajectory;Forces} is the probability of realiza- 
tion of the time-reversed process under time-inverted perturbations and conditions. The 
operation ~ means the inversion of time direction and besides of signs of variables which 
have the odd time parity. From here, after integration over all possible processes, the re- 
markable equality (exp(— Work/T)} = 1 does follow what leads to inequality (Work) > 0. 
Particularly, if the perturbation is described by Hamiltonian 

H(q,p,t) = H (q,p) - J2Qj( ( l^P) x j( t )^ 

with q, p and x(t) being canonical coordinates and momentums and external forces, 
respectively, one has Work = J £ Jj(t)xj(t)dt , where Jj(t) = f t Qj(q(t),p(t)) are the 
flows (velocities, currents, etc.) conjugated with the forces. 

This FDR is valid also in quantum case, under certain definition of quantum character- 
istic and probability functionals [34]. Besides, this FDR can be extended to nonequlibrium 
initial ensembles [34,38], in other words, to a combination of dynamical perturbations 
(i.e. thats of Hamiltonian) and thermal perturbations (i.e. thats of probability distri- 
bution in phase space). The producung FDR looks quite similar besides that the factor 
exp(— Work/T) should be replaced with exp(— AS) where AS'jt, Trajectory; Forces} has 
the sense of increment of entropy (or similar thermodynamic potential) during Trajectory. 



B. FLUCTUATIONS OF DISSIPATION 



Let ( J(t±), J(t n )) q denotes the n-th order nonequilibrium cumulant corresponding 
to the time-cut modification of perturbing forces, x(t) =>• x(t)r](t — min(ti, t n )) , where 
rj(t) — 1 at t > and rj(t) = at t < (thus the most early variable belongs to yet 
equilibrium state). Such correlators were termed in [34] quazi-equilibrium. Then the 
rigorous relation takes place [33,34,11], 



what extends Kubo formulas to arbitrary non-linear responce. As a consequence, the 
average value of the work (i.e. of the energy absorbed by the system), (Work) = 
f-oo (■AO) x (t)dt , can be simply expressed via two-point quazi-equilibrium correlators. 

The similar formulas for higher cumulants can be obtained [34] which lead to the 
formally exact expression for variance of the work [11] : 

(Work {2) ) = 2T(Work) + - J* x(l)x(2) (J(l), J(2), J(3)) q x(3)dld2d3 

where (Work^^ = (Work 2 ) — (Work) 2 , and for brevity the time variables (and indices) 
are denoted by integers and triple time integral is taken under condition 1 > 2 > 3 . 
Here the first term on the right-hand side describes shot noise accompanying the average 
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energy influx, and the second term involves an excess noise. The latter enters by means of 
three-point cumulants, so one could not snatch at it in Gaussian approximation. Under 
infinitesimally small perturbation, x — * 0, the excess contribution is of order of x 4 and 
can be expressed as 

(Work {2) ) = 2T(Work} + 



+- f x(l)x(2) (J(l), T(2, 3), J(4)) x(3)a;(4)dld2cZ3d4 + 0{x & ) 

1 J —oo 

where the dynamic linear differential responce function is introduced, 

T 3k (t,t')^[5J 3 (t)/5x k (t% =0 

the subscript denotes averaging over equilibrium Hibbsian ensemble, and integrations 
are ordered by the conditions 1 > 4, 2 > 3, 2 > 4, 3 > 4. 

If low-frequency fluctuations of kinetical values responsible for energy dissipation take 
place then these fluctuations are reflected just the excess noise term. In a steady state, 
when (Work) oc (t — t )x 2 , with to being the start time of perturbation, a l// 7 -noise 
should result in oc (t — t ) 1+1 x 4 behaviour of the excess part of energy variance. Hence, 
nearly equilibrium fluctuations of dissipative values (resistance, friction, mobility, etc.) 
can be generally reduced to four-point equilibrium cumulants like (J(l), T(2, 3), J(4)) 
which describe specifically non-Gaussian correlations between random responce and noisy 
flows. In their turn, these correlators can be connected with literally four-point cumulants 
(J(l), J(2), J(3), J(4)) , due to the four-point FDR [32,11], in part due to 

i(J(l),J(2),J(3),J(4)) = 



5x{A) 



(J(1),J(2),J(3)> 



x.=n 



= ( J(l), r(2, 4), J(3)) + (r(l, 4), J(2), J(3)) + (J(l), J(2), T(3, 4)> 

where the time moments are ordered as 1 > 2 > 3 > 4. More rich information is given by 
the 7 independent four-point relations derived in [32], in particular, by the FDR 



1 

2^2 



(J(l), J(2), J(3), J(4)) = 5 2 (J(l), J(2)> /Sx(3)5x(A) 



x=0 



+ 



+ 



5 2 (J(2), J(4)) /5x(3)5x(l 



x=0 



where the same time ordering is supposed, 1 > 2 > 3 > 4, and tilda means Xj(tj) = 



SjXj ( 



£j Jj ( 



-tj) , with Sj being the time parities of forces. 
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C. CONNECTIONS BETWEEN EQUILIBRIUM 
AND NON-EQUILIBRIUM RANDOM WALKS 



Let Q(t) = Q(q(t),p(t)) be coordinate of some random walk and x(t) be generalized 
force conjugated with Q in the above defined sense. In two of many possible examples 
Q(t) is a coordinate of Brownian particle with x(t) being the external force acting on it, 
or Q(t) is electric charge transported through some conductor with x(t) being applied 
voltage. Let Q(t) be time-even variable, the force remains constant after its switching on 
(x = 0att<0,x = const ^ at t > 0) and Q — at t — . Then the general 
producing FDR results in simple relation [7,9,38], 

W(Q,t)eM~^) =W(-Q,t) 



W(Q, t) - W(-Q, t) = [W(Q, t) + W(-Q, t)} tanh(^) 

where W(Q,t) is the probability density of Q at t > 0. From here one gets (Q(t)) = 
(Q(t)tanh[Q(t)x/2T]) . The expansion of the latter relation over x up to x 3 yields 

im)) * 



x=0 



(Q 2 (t)) /2T,and [9] 



dx 3 



(Q(t)) 



x=0 



3 

2T 



dx 2 



(Q(t),Q(t)) 



x=0 



AT 3 



(Q(t),Q(t),Q(t),Q(t)), 



with T being the initial temperature 

The first of these formulas is nothing but Einstein relation between low-field mobility 
and equilibrium diffusivity, or Nyquist relation, etc. The second formula asserts the 
connection between quadratic non-equilibrium addings to variance of Q(t) (that is the 
excess noise), equilibrium fourth cumulant and cubic response. The latter could be due to 
either nonlinear dissipation mechanisms or to dependence on actual current temperature 
maintained by Joulean heating. But, in any case, the left-hand side can not grow in a 
more fast way than oc t. On the contrary, in presence of /~ 7 excess noise the first term on 
right-hand side should grow as oc t 1+7 . Hence, at sufficiently long time intervals it must 
be compensated by the second term, and we come to relations between low-field excess 
nonequilibrium noise and four-point cumulant of equilibrium noise. The result can be 
written also as K m {r) = jvK d (t) [4,6,9], where K d (t) is the above defined correlator of 
diffusivity and K m (r) is the mobility correlation function as introducedd by the expansion 
(Q(t), Q(t)) = (Q(t), Q(t)) + x 2 ft / * Kjt' - t")dt'dt" + 0(x 4 ). This is the extension of 
Einstein relation and similar relations to 1/ f-fluctuations. 



D. UNIFICATION OF EQUILIBRIUM AND 
NONEQUILIBRIUM BROWNIAN MOTION 



The producing FDR can be reformulated in terms of characteristic function (CF): 



x 



exp[(ifc - -)Q(t)] = (exp[-tkQ(t)}) 
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This is functional equation with respect to dependencies on ik and x . It helps to find 
the connections between the whole statistics of nonequilibrium random walk and that of 
equilibrium one, in presence of 1/f- noise, if we know something about statistics under 
neglecting 1/f- noise [4,6,8,38]. 

One of possible solutions on the functional equation is (exp[ikQ(t)]) = E(ikx/T — k 2 , t) 
, with only two independent arguments. It is easy to see that both the CF of quazi- 
Gaussian Brownian motion and CF of asymmetric one-directed random walk with Cauchy 
statistics can be unified into just such form [6,8], 

<exp{ifcAQ}> *Z(ikx/T-k 2 ,t) , E(£,t) = exp{-D ^tA(t)hi(j ~ ^)) 

Here x = corresponds to quazi- Gaussian case, while the change ikx/T — k 2 =>- ikx/T 
, under the condition \k\x/T » k 2 , leads to the second case. The condition means 
that only the drift component of the walk is observed which dominates at large time 
scales. Comparing the CFs we find that the unification needs in relations U = D x/T 
and Q = q^x/T (or V = D'x/T ). These is nothing but Einstein formula (or Nyquist 
formula, etc.) and its generalization for 1/f-fluctuations. 

Besides, we obtain the statistical model of nonequilibrium Brownian motion, with tak- 
ing into account both the drift and diffusional displacements and common 1 / f-fluctuations 
of diffusuvity and mobility (in [6] FDR were used even in order to derive the non- 
equilibrium Cauchy statistics from quazi- Gaussian statistics obtained in [4,5]). The ad- 
ditional generalization of the model allows for arbitrary dependencies D = D (x) and 

jy = d'(x). 



VII. CURRENT STATE AND PARADOXES 
OF GAS KINETICS 

A. BBGKY HIERARCHY 

The probabilistic formulation of the Hamiltonian dynamics of classical gas is given by 
Bogolyubov-Born-Green-Kirkwood-Yvon (BBGKY) equations 



d_ d_ 

dt Vl dr\ 



F 1 = p J A 12 F 2 dp2dr 2 = I 12 (F 2 



d d d . 

at or i or 2 



F 2 = pj (Ai 3 + A 23 )F 3 dp 3 dr 3 = I 13 (F 3 ) + I 23 (F 3 ) 



d_ 

dt 



F n — P J Kj n +iF n+ \dp n+ idr n+1 = y~^ j J,- ra+1 (F Tt+1 ) 



dr 3 , ; 
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Here is the n-particle Liouville operator, F n = F n (t,ri..r n ,pi..p n ) is the n-particle 
probability distribution function (DF), normalized to whole (formally infinite) gas volume 
Q, J F n dp n dr n /VL = F n _i, the operators Ajk decribe pair interactions, with f(r) being the 
interaction force, and Vj = Pj/m with m being the mass of particles. 

B. BOLTZMANN EQUATION 

More than 100 years ago Boltzmann suggested the phenomenological kinetical equation 
to describe the evolution of rarefied gas with a short-range repulsive interactions between 
particles. The Boltzmann equation (BE) can be written as 



d d 
di + Vl dR 



= pJd P2 J d 2 B \ Vl - v 2 \ {Ft{t,R,^) - F* n (t,R,p u p 2 )} = S 12 (F?) (1) 
where the ansatz about molecular chaos (the famous Stosshalansatz) is introduced, 
F* n {t,R, Pl ,p 2 ) = F l {t,r l w i2,pi)Fi(i,r 2 « R,p 2 ) 

Here F™(t, R,pi,p 2 ) is DF of two particles which come into collision at a point R from 
the preceeding in-state, B is the two-dimensional impact parameter vector perpendicular 
to the relative velocity v — v\ — v 2 ( B = r — v(v , r) / \v\ 2 with r = r\ — r 2 ), and p[,P2 are 
the momentums what should be attributed to initial in-state in order to transform into 
the momentums pi,p 2 at final out-state after collision. In accordance with Stosshalansatz, 
the join two-particle DF is factored into the product of one-particle DFs. 



C. IS BOLTZMANN EQUATION DERIVABLE 
FROM BBGKY HIERARCHY ? 

With no doubts, BE is the well grounded kinetical model. Thus it is the striking 
surprise that it still is not derived from BBGKY equations. 

The widely spread opinion is that the derivation was performed by Bogolyubov. How- 
ever, it was based on two ansatzes. First, the evolution of statistical ensemble governed by 
BBGKY equations tends to the stage when all the DFs become time-local functionals of 
one-particle DF. Secondly, the molecular chaos (mistakenly identified wih the weakening 
of correlations of far spaced particles). 

Both these assumptions seem quite reasonable. But if it is really so, both should be 
deducable from BBGKY formalism. The possibility of rigirous derivation of BE from 
BBGKY equations was carefully analysed by O.Lanford, under a suitable definition of 
the norm in DF's functional space in the Boltzmann-Grad limit [35,36] (see also [30,37]). 
It was shown that the divergency of the series expansion of solution of BBGKY hierarchy 
to solution of BE takes place, but only during a short time after start smaller than the 
free flight time. 
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The attempts to extend this result to larger time intervals run into likely principal diffi- 
culties. The matter is that the singularities (discontinuities) of binary DF and higher DFs 
arise at hypersurfaces in the phase space which correspond to collisional configurations 
(for instance, when the vector v 2 — V\ is parallel to r 2 — r\ ). Though these singularities are 
concentrated in regions whose Lebesgue measure tends to zero in Boltzmann-Grad limit, 
thats are just the regions which determine the evolution of DFs of colliding particles and 
reliability of molecular chaos. 

D. OBJECTIONS TO MOLECULAR CHAOS 
AND BOLTZMANN EQUATION 

M.Kac [31] had doubts that BE is derivable from the BBGKY theory in general 
spatially non-uniform case. He treated BE sooner as the good model than as formal 
approximation of exact theory. According to Kac, in the uniform situation one could 
come to BE by means of spatial averaging in addition to ensemble averaging. In this 
case DFs contain no information about previous motion of gas particles, therefore the 
statistical ensemble makes no difference between two actually colliding particles and two 
arbitrary particles, and molecular chaos looks quite inevitable. But in general one losses 
the help of spatial averaging. 

This criticizm was developed in [9]. According to BE, a life of a gas particle is pre- 
scribed by quite certain portions of collision of various sorts (by the portioning what 
corresponds to uniform distribution of the impact vector B ). But in reality the distribu- 
tion over sorts is out of both dynamical and thermodynamical control of the system. Even 
if it happened that some portion was much larger than its expected value, this incident 
will not influence the future. There are no mechanism what would be able to remember 
the history of preceding collisions and use it in order to compensate previous occasions 
and keep a certain distribution (histogram) of collisions over sorts. 

Consequently, in accordance with the above phenomenology, the number of collisions 
(of any given sort) per unit time may have no predictable limit. In other words, the rate 
and the efficiency (the effective cross-section) of the flow of collisions of a gas particle 
undergo scaleless flicker fluctuations. The effect of these fluctuations on the evolution of 
Fi could be diminished only by spatial averaging in the uniform gas, but in no way in 
the framework of non-uniform ensemble, when the spatial dependencies of DFs carry an 
information about history of Brownian motion of gas particles [9]. 

Indeed, the rate (the diffusivity) of the Brownian motion is closely correlated with 
fluctuations in the number and efficiency of collisions, hence, with the currently possible 
collision too. Therefore, the present collision has long-living statistical correlations with 
the past motion of colliding particles. In view of these correlations, the binary DF for 
particles going into collision acquires the specific statistical meanig as the conditional 
probability distribution, under the condition that the joint collision is realized in a fixed 
space-time point. Hence, this DF can not be factored into product of DFs F\. 

The generally possible factoration has the form 

Fi n (t,R, Pl ,p 2 ) = F{(t,R, Pl )F{'(t,R,p 2 ) 
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with conditional marginal distributions different from the conditionless one-particle DF 
Fi, F[ ^ Fi, F" 7^ Fi , what corresponds to mutual independency of particles in the 
momentum space but not in the real configurational space. Thus, the molecular chaos 
fails, at least because of spatial statistical correlation of colliding particles. 

Such the ensemble can not be reduced to BE. It is amusing that the violation of molec- 
ular chaos takes place only between colliding particles, that is just where the molecular 
chaos is especially believed. This picture is obviously incompatible also with the ansatzes 
used by Bogolyubov. 

Therefore, we must introduce a different ansatz, with the same purpose to build a ki- 
netical model and firstly to narrow a variety of kinetically describable solutions of BBGKY 
equations. 

VIII. CORRECT REFORMULATION 
OF BOLTZMANNIAN KINETICS 

A. PRIMARY ANSATZ OF GAS KINETICS 

The essence of coarsened kinetical language is that all the set of successive dynamical 
stages of a two-particle interaction process is treated as the unique momentary collision 
event. This event looks as a "black box" (compressed into point under Boltzmann-Grad 
limit) which instantly transforms the input velocities into output ones. 

Clearly, this collision event should conserve probabilities, else a model would effectively 
allow for creation or annihilation of particles inside the black box. In other words, in any 
kinetical statistical ensemble, the flow of particles entering the collision should currently 
coinside with the flow of particles what leave it. 

In our opinion, just this trivial requirement should serve as the basic ansatz when de- 
riving a kinetical model from BBGKY equations. In order to provide with the probability 
conservation, the exclusion of internal dynamics of collision from consideration must be 
compensated by definite bondary conditions at the black box borders. Concretely, one 
should keep the probability measure of any final output state at the border to be currently 
equal to the probability of corresponding initial input state. 

This condition is formally similar to the one for DFs at \rj — Tk\ — d in the rigid sphere 
gas (with d being the diameter of spheres). In general, a choise of the borders of collision 
box is rather arbitrary, at least in the framework of Boltzmann-Grad limit. Therefore 
the above mentioned ansatz is equivalent to that probabilities of all the inner stages of 
any fixed sort of collision equal to one and the same quantity. If this condition was not 
satisfied, it would be impossible to use the concept of momentary collision event at all! 

In such the way we narrow the set of solutions of BBGKY hierarchy to be under 
kinetical consideration: the distributions somehow modulated along a path inside collision 
box are excluded (although arbitrary modulation from one path to another is permissible). 

The full characterization of a sort of binary collision consists of its location defined by 
three space coordinates, the two-dimensional impact parameter vector and two velocities 
of particles, that is totally by 3+2+6=11 scalars instead of total 12 phase space variables of 
the pair. The rest twelfth variable, let be denoted as , just enumerates inner dynamical 
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stages of the collision in the black box interior. This variable should be excluded from a 
kinetical picture. Thus, with respect to F 2 , the primary ansatz means that 

d 

— -F 2 = , r 2 — r\ G collisionbox (2) 

One always can choose 6 as the inner time of the collision, for instance, as the time 
necessary to achieve its given stage from its most proximity stage. Then the differentiation 
by looks as 

with r = r 2 — r\ and v = v 2 — v i being the relative distance and velocity of particles, re- 

(2) 

spectively. The operator L T ^ is nothing but the part of full two-particle Liouville operator 
l/ 2 ) responsible for the inner dynamics of collision. 

Of course, none of rigorous solutions on BBGKY equations exactly satisfies (2), as 
well as the ansatzes introduced by Boltzmann and Bogolyubov. But if the kinetic ap- 
proximation is actually available, we may hope that at a more or less late stage (kinetical 
stage) of the ensemble evolution the ansatz (2) is quantitatively satisfied to the extent 
sufficient to regard it as formal identity. 

B. COLLISION INTEGRAL 

In fact, one inevitably needs in the primary ansatz in order to construct the Boltz- 
mannian collision integral, even under the molecular chaos postulate, and thus to reduce 
the explicite dynamics to the momentary collisions representation. Indeed, let us rewrite 
the right-hand side of first BBGKY equation as 

^12(^2) = P J A 12 F 2 dp2dr 2 = 

= P S L {V2 - V ^r F2drdp2 - P J Lb Ij§ F ' drdp2 

where CB marks the spatial integration over the collision black box. After transformation 
the volume integral in the first right-hand term into surface integral § over the collision 
box board and after separation of in-states and out-states, regarding of sign of the scalar 
product (vr) , one obtains 

h2(F 2 )=pJ 

Here in both the surface integrals one and the same velocity arguments are present while 
the impact parameter vectors B' and B relate to final and initial state, respectively. 

This expression can be reduced to conventional collision integral only if the probabil- 
ities of the final states constantly coinside with thats of preceding incoming states, 



pout d 2 B f 



vr>0 
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F™%R, Pl ,p 2 ,B') = F™(t,R,p[,p' 2 ,B) 

Here R is the location of collision, and in view of Boltzmann-Grad limit we do not make 
difference between spatial locations of in-states and out-states. 

The validity of this relation is ensured just by ansatz (2). Besides, due to (2) the last 
integral, what formally acts like a bulk source of particles, turns into zero. As the result, 
we come to the kinetic equation (1), but, as it will be seen, the molecular chaos hypothesis 
fails. 



C. VIOLATION OF MOLECULAR CHAOS 

According to primary ansatz, in the interior of collision box the second BBGKY equa- 
tion reduces to 



d d d 

di + 2 dR + dO 



d 8 

m + U2 dR 



F 2 = I 13 (F 3 ) + I 23 (F 3 ) (4) 



with R = (ri + r 2 )/2 and U 2 = (t>i + v 2 )/2 being the position and velocity of the mass 
center of colliding pair, i.e. the position and velocity of the collision as a whole. 

Hence, the spatial drift of the DF of colliding or closely spaced particles is governed 
by the mass center velocity. This consequence of (2) has very simple physical meaning. 
First, if the relative motion of particles is included into the collision interior, then quite 
naturally it should be automatically excluded from the transport of the collision as a 
whole. Secondly, the DF related to collision box represents the conditional probability 
under the condition that the particles are partners in one and the same encounter. It 
is natural that this specified DF is ruled by an equation different from what describes 
the less informative unconditional DF. Evidently, the conditional DF F 2 in (4) can be 
interpreted as the probability measure of that a collision occures in the point R , i.e. as 
the measure of ensemble-averaged density of collisions. 

It is easy to see that in a spatially non-uniform ensemble, when 8F 2 /dR ^ , the Eq.4 
forbids the factoration of this DF into product of one-particles DFs, F 2 ^ F\(t,pi, R) * 
Fi(t,p 2 , R) , regardless of contribution of the right-hand side. Therefore the molecular 
chaos inevitably fails, and just with respect to colliding particles! As a consequence, in 
general the density of collisions is not simply proportional to [density of particles] 2 . 



D. 3-PARTICLE PROCESSES AND EVOLUTION 
OF DENSITY OF COLLISIONS 

Let us take in mind the Boltzmann-Grad limit (BGL), 

p — > oo, d — > 0, A = 1/pd 2 = const 

where A is the mean free path and d is characteristic radius of interaction. 

The right-hand side of the Eq.4 describes a change of density of collisions due to three- 
particle processes. Thats are not factual triple collisions, which are of no importance under 
BGL, but two successive pair collisions separated by a distance comparable with d but 
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perhaply much greater than d . One of the collisions involves a third external particle 
representing the rest gas. Under BGL, a spatial region taken by this process contracts 
into point if measured by A scale although can grow up to infinity in d units. 

In order to express the right-hand side of Eq.4 in kinetic fashion, in terms of collision 
integrals, we need in the obvious three-partricle extention of the ansatz (2), 

dF 3 /dO = L®F 3 , L?J = - U 3 d/dR 

where 6 is a suitably choosen inner time of the process, and R and U 3 are the position 
and velocity of mass center of the close three-particle configuration under consideration. 
Eventually, the Eq.4 is transformed into the kinetical equation 

F 2 = S 13 (F™) + S 23 (Ff) (5) 

where the superscript in means that the external third particle is at in-state with respect 
to 1-st and 2-nd particles (i.e. is coming to them from the infinity, in d units). This 
equation describes the evolution of density of binary encounters and collisions as changed 
by collisions of the given two particles with the rest gas. 

E. NO FAR CORRELATIONS 
ALTHOUGH NO MOLECULAR CHAOS 

In contrary to the conditional binary DF subjected to Eq.4, the general binary DF for 
arbitrary configurations placed outside the collision box is subjected to the usual equation 

F 2 = S 13 (FD + S 23 (FD (6) 

Evidently, this equation allows for the factored solutions, F 2 = F 1 (t,r 1 ,p 1 ) * F 1 (t,r 2 ,p 2 ) 
, with F™ = F™(t,r 1 ,p 1 ,p 3 ,B)*F 1 (t,r 2 ,p 2 ) in S 13 (Ft) and F 3 m = F™(t,r 2 ,p 2 ,p 3 , B) * 
Fi(t,ri,pi) in S 23 (F 3 n ) , where each of F\ multipliers is the solution on the Eq.l. 

Thus the vanishing of correlations of far distanced particles coexists with invalidity of 
molecular chaos with respect to closely spaced particles (in particular, colliding ones). If 
measure the inter-particle distance with the free flight path A then, under BGL, the weak- 
ening of correlations, F 2 = Fi(t,ri,pi) *F\{t, r 2 ,p 2 ) , can be rigorously satisfied anywhere 
except the hypersurface r 2 — r\ = . In view of this circumstance, the singularities of 
DFs at collisional hypersurfaces met in [35,36] confirm the necessity to specially consider 
the specific DFs for encountering configurations with r\ = .. = r n . 

F. MANY-PARTICLE ENCOUNTERS AND 
CHAIN OF KINETICAL EQUATIONS 

In the framework of BGL, one should continue the above mentioned procedure up to 
infinity, with extention of primary ansatz to many-particle events, namely, to sequences 
of pair collisions and close many-particle encounters, step by step building the kinetical 
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model of spatially non-uniform gas. The extention looks as dF n /dQ = , with the 
equality to be satisfied inside a collision (encounter) black box alloted for the event and 
being its inner time. This is the inevitable payment for decision to reformulate the 
ensemble evolution in terms of collision integrals. 

Because of identities d/dQ = L^j , = U n d/dR + L^j , where L^j is the part of 
n-particle Liouville operator responsible for relative motion inside the event, the result is 
the infinite chain of coupled kinetic equations [9] 

n 

i^Ew^r+i) (7) 

Here R = (r ± + .. + r n )/n , U n = (vi + .. + v n )/n . 

The DFs on the left represent the probability measure of occurence of n-particle en- 
counters, in particular, constituted by (n — 1) successive collisions. The right-hand DFs 
relate to (n + l)-particles events, with an external (n + l)-th particle being at precolli- 
sional in-state with respect to the left n-particle configuration. In general, any of these 
F n depends on time, the location of the encounter R , n momentums, (n — 1) impact 
parameter vectors and, besides, on (n — 2) distancies between binary encounters if n > 2 . 
Again, the relative motion turns into the hidden inner property of the event and so does 
not influence the drift of the probability of the event as a whole. 

In spatially uniform case, when d/dR — + , the solution on the Eqs.7 can be sought in 
the conventional molecular chaos form, F n = Y[j F±(t,rj,pj) , with no factual dependence 
on impact parameters, and the model reduces to conventional BE. 

G. MORE ABOUT PRIMARY ANSATZ 

To formulate the primary ansatz specially for BGL, let us write the interaction force 
in the form /(r) f(r/d)/d and measure the position of a cluster of close particles in A 
units (and time in free flight time units) while interparticle distancies in d units. Then 
BBGKY equations turn into 

= / Y,^n + iF { n %d Pn+l dr n+1 
J j 

where /i — d/X — pd 3 — > is small gas parameter, and DFs dependence on /i is marked. 
Clearly, if a definite limit of all the does exist at pi — > , then without fail the relation 
dF^/dQ = L^F^ = is satisfied. Hence, the ansatz is equivalent to the assumption 
that L^F^) = o(/i) , at least for configurations corresponding to chains of far spaced 
binary collisions (i.e excluding actually triple or multiple collisions). 

H. COLLISIONS, ENCOUNTERS AND 
COARSENED KINETIC EQUATIONS 

In principle, under a suitable geometrical parametrization of the events the left-hand 
DFs in Eqs.7 can be identified with the right-hand ones. However, a full classification of 
configurations would lead to an extremely complicated theory. 
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It seems reasonable to simplify the model with the help of spatial averaging of DFs 
over Tj inside the collision (encounter) regions, under fixed mass center position R . Of 
course, then a dependence on impact parameters will be lost, except the only parameter 
on the right-hand sides what relates to the incoming (n + l)-th external particle. As a 
consequence, one will be forced to distinguish left-hand DFs F n and right-hand DFs F™ . 

Besides, for more simplicity, in the framework of BGL we may choose the dimensions of 
the collision (encounter) regions (boxes) much larger than interaction diameter d . Then 
the most part of configurations under averaging will correspond to encounters, when 
particles are separated by a distance comparable with d but without factual interaction 
(and so out-states coinside with in-states). 

We may hope that in spite of such the crucial simplification the essential features 
of statistics of collisions will not be lost. Indeed, in view of incident origin of close n- 
particle clusters, there are no principal difference between efficient collisions and merely 
encounters. The latters can be regarded as collisions with negligibly weak scattering, and 
fluctuations in the flow of collisions should be similar to thats in the flow of encounters. 

I. WEAKENED MOLECULAR CHAOS 

It is interesting that the chain (7) is still formally time-reversible if all the geometry of 
binary and many-particle events is kept under control. Indeed, the combined inversion of 
time and velocities and rearrangement of in-states and out-states does not change these 
equations. However, this does not prevent existence of seemingly irreversible solutions 
of the equations, all the more if a part of the complete characterization of many-particle 
configurations is lost. 

The DFs smoothed over collision regions have the clear physical meaning: thats rep- 
resent ensemble-averaged densities of n-particle encounters, each depending on time, lo- 
cation R and n velocities. But in order to close the Eqs.7 one needs in some connections 
between F%] rl and the left-side functions. This is just the point where one more ansatz 
like molecular chaos should be attracted. 

It is natural to assume that the velocity distribution of (n + l)-th right-hand external 
particle is statistically independent on velocities of the left-hand cluster particles: 

f£.iM>P (n+1) ) = G B (f,i2,p< B >)F 1 (f,i2 J p B+1 ) 

where the notation p^ = {p±, ..,p n } is introduced. In general this factoration is not iden- 
tical to absolute statistical independence of particles, since it allows for spatial correlations 
whose manifestation is that functions G n differ from F n . 

Clearly, one should connect F%] x and thus G n with F n+1 , not with F n . The matter 
is that the spatial dependencies of both F™_ x and F n+i reflect the probability density of 
the same (n + l)-particle encounter events. Therefore the degree of specifically spatial 
correlations must be one and the same in both F%]_ x and F n+1 (this assertion follows also 
merely from the proximity of (n + l)-th particle to the left-side cluster in A units). This 
requirement implies 

J F n+1 dp( n +V = J F™ +1 dp( n+ V J G n (t,R,pW)dpW J F 1 (t,R,p')dp' 
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One can say also that this relation expresses the conservation of particles and probabilities 
during encounters. 

Besides, take into account that the degree of correlations inside the left-hand cluster 
must be the same in F^" x and in F n+1 . Then the only reasonable form of G n is G n = 
J F n+ idp n+ i/ J Fidpi, and we finally come to the connection [9] 

This is the weakened form of molecular chaos which establishes the statistical inde- 
pendence of the external particle only in the velocity space, but not in the configurational 
space. 

J. NON-UNIFORM EQUILIBRIUM GAS 
AND BOLTZMANN-LORENTZ MODEL 

One naturally is forced to consider some spatially non-uniform solutions on BBGKY 
hierachy even in thermodynamically equilibrium case, if wants to analyse the statistics 
of Brownian motion of a test gas particle. If the position of even a single test particle 
is known at initial time moment, or even merely its distribution differs from the uniform 
one, then one should deal with spatially varying DFs. 

Let DF Fi(t,R,pi) be related to the test particle. Then the higher order DFs F n+1 
subjected to Eqs.7 and 8 relate to the clusters consisting of the test particle plus n 
arbitrary close particles. All the DFs together encode the information about statistics of 
encounters, collisions and random displacements of test particle. 

But, of course, now the DF of test particle, F\(t,R,pi) , and the one-particle DF 
Fi(t, R,p n+ i) in (8) are different things. The latter relates to arbitrary gas particle and 
so should coinside with the Maxwellian distribution in equilibrium gas: 

F 1 (t,R,p n+1 ) = F e «(p n+1 ) , F e «(p) = (27rTm)- 3 / 2 exp(V/2mT) 

The weakened molecular chaos condition transforms into 

F™ +1 (t,R,pW,p n+1 ) = F e «(p n+1 ) J F n+l (t,R,p( n \p>)dp' (9) 

The sets of equations (7) and (9) represent the analogy of the well known Boltzmann- 
Lorentz model but correctly modified for a spatially non-uniform ensemble. 

IX. SELF-DIFFUSION AND 
DIFFUSIVITY 1/F-NOISE 

A. KINETICAL EQUATIONS FOR BROWNIAN 
MOTION OF TEST PARTICLE 

Let r(t) and u(t) be coordinate and velocity of any particle and Q(t) = ^u(t')dt' 
be its displacement during time interval t (more accurately, let u and Q(t) be the pro- 
jections of velocity and displacement vector onto a fixed direction). The statistics of 
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self-diffusion (Brownian motion) of gas particle is described by the characteristic function 
(CF), (exp[ikQ(t)]) , but, as we know, the most significant information is contained in the 
first two terms of expansion of In (exp[ikQ(t)]) over k 2 , which reflect the mean diffusivity 
and the correlation function of diffusivity fluctuations. 

The Laplace transform of CF always can be represented as 

roo 

/ (exp[ikQ(t)]) exp(-zt)dt = [z + kD(z, ik)k]~ l 
Jo 

with D(z, ik) being the frequency dependent diffusivity kernel. From the other hand, with 
the help of generalized Green-Kubo formulas, one can connect D(z,ik) and the response 
of probability distribution of a test particle to a weak external force / ex t(r) , f ex t , 
which influence it after switching on at t = (being parallel to the choosen direction). 
Concretely [9], 

roc r r 

J dtexp(-zt) J drexp(-ikr) J vF 1 (t,r,p)dp — 



D{z,ik) 



T[z + kD(z,ik)k] 



exp(-ikr)f ext (r)dr 



where v — p/m , the DF Fi(t,r,p) describes the test particle and equals to F eq {p) at 
t < . 

Hence, the long-time statistics of equilibrium self-diffusion can be obtained from the 
spatially non-uniform evolution (formally nonequilibrium but in fact equilibrium in the 
macroscopic sense) of ensemble disturbed by a force what acts on the only test particle. 
In the framework of the kinetical theory constructed, the problem reduces to solution of 
Eqs.7 and 9 , after adding the terms f ext (R)dF n /dpi to all left sides of Eqs.7. 

We may avoid a part of non-principal mathematical difficulties if replace the Boltzman- 
nian collision integrals on right-hand sides of Eqs.7 with the Fokker-Planck approximation 
of the Boltzmann-Lorentz collision operator, 



S jn+1 ( I F n+1 (t, R,p^,p')dp' * F^{p n+l )) => 7 — ( Pj + mT—) j F n+1 dp n+l 



with 7 being the inverse free flight time. Then we have to solve the equations [9] 

n d d r 

F n = Yl 'y-g—. (Pj + mT ~Q~) J F n+l d Pn+l (10) 



d d d 

di + Un dR + fext{R) dp~i 



i=i d Pj d Pj 
with U n = (vi + .. + v n )/n and equilibrium initial conditions F n \ t =o— IT; F eq {.Vj) ■ 

B. NON-LOCAL DIFFUSIVITY KERNEL 

But even the analysis of these simplific equations is rather hard, and only the lowest 
terms of the expansion 

D(z,ik) = D (z) + {ik) 2 D 1 {z) + ... 
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were found in [9] . Nevertheless, these terms always allow to check the existence of low- 
frequency duffusivity fluctuations and estimate their variance. 

In accordance with the D(z,ik) definition, one can verify that D (z) is the Laplace 
transform of the two-point equilibrium velocity correlator, K u (t) = (u(t)u(0)) , and 

rco rt 

D 1 (z) = / {K D (t) + K u (t) / (t-t')K u (t')dt'}exp(-zt)dt 
Jo Jo 

Here Ko{t) is the above discussed diffusivity correlation function as expressed via four- 
point equilibrium cumulant of u(t) . Therefore, at low frequencies, \z\ << 7 , D\{z) is 
nothing but the Laplace transform of K D (t) , that is Re[Di(2irif + 0)] is responsible for 
the power spectrum of diffusivity fluctuations. From the other hand, one can connect 
Di(z) strictly with the fourth statistical moment of the displacement: 

roo , . 24 

j Q (Q\t))eM-zt)dt =-[ Dl (z)+D 2 (z)/z} 

The latter function was obtained in [9] in the form of double series via 7/2 which can 
be exactly transformed into an integral. Here we omit these manipulations. At \z\ << 7 
, the final result reeds as 

1 t 

D i( z ) ~ ir D o( z )\^ 2 - + const \ 

2z z 

with D (z) = ^(7 + z)~ l and D = D o (0) = T/mj being the usual diffusivity. In the 
temporal reprezentation, these results mean that (Q 2 (t)) = 2D t and 

(Q (4) (t)) = (Q\t)) - 3 (Q 2 (t)) 2 » 6L> 2 t 2 [ln 2 7 t + c'ln 7 t + c"} 

at yt » 1 , where c', c" are numerical constants. In view of the above consideration, this 
is the clear evidence of existence of 1/f type self-diffusivity fluctuations. 

Formally, these results are produced by the peculiar drift terms U n -J=^ on right-hand 
sides of the Eqs.10. But our experience showed the remarkable fact that the concrete 
form U n — (vi + .. + v n )/n is of no principal importance. Almost any sequence of U n = 
CnjVj with different c n j results in an excess contribution to the fourth cumulant (^Q^ (t)^ 
which grows faster than proportionally to time and thus corresponds to a low-frequency 
diffusivity noise. In particular, if the mass of the test particle M differs from masses 
of gas particles, and so U n = (Mv\ + mv 2 + ...)/(M + nm) , then the Eqs.10 lead to 
(Q^\t)) ~ t 1+1 > 7 > , where 7 differs from unit and depends on the ratio M/m [9]. 

C. 1/F-FLUCTUATIONS OF DIFFUSIVITY 

The power spectrum of the relative diffusivity fluctuations resulting from the obtained 
Di(z) for sufficiently low frequences, / << 7 , is 
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Formally this is non-integrable spectrum like that of a slightly non-stationary noise. But 
we already emphasized that it is fictituous non-stationarity what refers only to the dura- 
tion of observation, not to the time moment when observation takes beginning. Besides, 
it is necessary to underline that the spectrum is insensitive to spatial dimensionality of 
gas, giving the illustration of so-called zero dimensionality of 1/f-noise. 

X. RESUME 

We reviewed the explanation of 1/f-noise based on the idea that its sources are not some 
specifically slow processes but, in opposite, the same processes which realize forgetting 
the past, relaxation and related white noise. If a system producing randomness constantly 
forgets its past and keeps in its memory scope only a limited amount of recently happened 
events, then it is unable to definitely distinguish a proper part of events and deviation 
from it. As a consequence, the fluctuations in the amount of incidently emerging events 
grow similarly to its average value, resulting in just 1/f-fluctuations of the rate of random 
flow of events. 

In order to comprehend such a random behaviour, one should make difference between 
long-living statistical correlations and long-living causal correlations. 

We noted that 1 / f-noises in different systems can be interpreted in such manner. The 
corresponding phenomenological models of purely random flow of events and of Brownian 
motion were considered, whose important feature is violation of the central limit theorem 
with respect to long-range statistics and close connections with Cauchy statistics. 

The attempts to develope an adequate philosophy of 1/f-noise were made. The quan- 
titative estimates of 1/f-noise in terms of only short-range (microscopic) characteristic 
scales were presented, showing that the wider is the memory scope of the system the 
lower is 1/f-noise level. 

For thermodynamical systems, the general connections between 1/f-noise and four- 
point equilibrium cumulants were analysed. The existence of long-living four-point corre- 
lations equivalent to flicker noise in exactly solvable Kac's ring model was found. 

The connections with Hamiltonian models of relaxation and irreversible processes were 
discussed, with hopes that if any realistic time-reversible many-particles dynamics pro- 
duces noise and irreversibility then it also produces 1/f-fluctuations of their rates. This 
expectation can be testified in such principally important cases as slightly unharmonic 
crystal lattice (phononic system) and classical gas with short repulsive interactions (in 
the Boltzmann-Grad limit). 

We pointed out that collisions of any gas particle with other particles form just such 
a random flow of events whose rate is being constantly forgotten and free of a thermo- 
dynamical control by the system. We argued that consequently 1/f fluctuations of the 
rate (number of collisions per unit time) and of related quantities (effective cross-section, 
self-diffusuvity, mobility) should occur and manifest themselves as the spatial statistical 
correlations of colliding particles in a spatially non-uniform statistical ensemble, resulting 
in violation of both the molecular chaos ansatz and the ansatz used by Bogolyubov to 
deduce Boltzmann equation from BBGKY hierarchy. 

The euristic reasonings were confirmed by critical analysis of kinetical theory of gas. 
It was argued that reformulation of the gas ensemble evolution in terms of coarsened 
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description with Boltzmannian collision integrals inevitably needs in another ansatz whose 
meaning is the conservation of particles and constancy of probabilities on the path what 
connects in-state and out-state through the collision box (i.e. the space-time region to be 
excluded from kinetical description). 

We showed that this requirement leads to the infinite set of specific kinetical equa- 
tions for the probability distributions on the hypersurfaces which correspond to binary 
and many-particle encounters and collisions (the sequences of binary collisions, not actu- 
ally multiple ones). In a spatially non-uiform ensemble these kinetical equations forbid 
conventional molecular chaos, because of producing statistical correlations of close parti- 
cles in the configuration space, and can be reduced to standard Boltzmann equation only 
in the limit uniform case. 

In our opinion, this is in agreement with the results by O.Lanford who noted the 
appearance of singularities of distribution functions at the above mentioned hypersurfaces 
in the formal series solution of BBGKY equations when investigating the possibility of 
reduction of BBGKY hierarchy to Boltzmann equation. The singularities just indicate 
that one needs in special consideration of probabilities of colliding and close configurations. 

The obtained chain of kinetical equations, with addition Green-Kubo formulas and 
higher-order fluctuation-dissipation relations, was applied to analyse the statistics of self- 
diffusion (Brownian motion of a test particle) in equilibrium gas. The four-order cumulant 
of Brownian displacement was calculated whose time behaviour proves the existence of 
1/f-fluctuations of diffusivities (and thus of mobilities) of gas particles. 

Hence, we have at least one example of principally correct formulation of kinetical 
approximation which does not loss 1/f-noise, and the example of dynamical system by no 
means producing slow processes but in spite of this producing 1/f-noise. 
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